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STATEMENT  OF  THE  PROBLEM 

The  phenomena  of  cavitation,  produced  in  liquids  by  an 
ultrasonic  beam,  was  observed  nearly  a quarter  of  a century 
ago.  It  was  noticed  that  many  previously  unobserved  effects 
accompanied  the  phenomena;  the  efficiency  and  transmitting 
power  of  the  transducer  as  well  as  the  physical  and  chemical 
properties  of  the  transmitting  medium  were  radically  altered. 
Many  of  the  other  physical  and  chemical  phenomena  produced  by 
ultrasonic  waves  were  found  in  some  way  to  be  connected  with 
cavitation,  £‘he  existing  theories  of  the  basic  structure  of 
liquids  were  soon  found  to  yield  information  which  was  not 
correlated  by  experimental  data. 

The  purpose  of  this  paper  is  to  make  a comprehensive 
survey  of  the  phenomena  of  ultrasonic  cavitation  in  liquids. 
For  better  understanding  of  the  problem  a brief  discussion  of 
the  general  properties  and  characteristics  of  ultrasonic  waves, 
as  well  as  a short  history  of  the  development  of  the  science 
of  ultrasonics,  has  been  included.  The  more  recent  methods 
of  production  of  high  frequency  sound  waves,  along  with  the 
experimental  methods  and  equipment  employed  in  the  study  of 
cavitation  has  also  been  briefly  described. 

It  is  the  main  purpose  of  this  study  (1)  to  review  the 
present  theories  of  the  ultrasonic  production  of  cavitation 
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in  liquids;  (2)  to  evaluate  the  theories  in  light  of  the 
experimental  methods  and  results;  (3)  to  analyze  the  prop- 
erties of  liquids  in  an  ultrasonic  field;  (4)  to  summarize 
the  more  tenable  theories  of  cavitation  in  view  of  present 
day  researches  in  the  field  of  ultrasonics. 

Although  it  has  been  impossible  to  thoroughly  review 
the  voluminous  literature  written  about  ultrasonic  phenom- 
ena, an  attempt  has  been  made  to  include  those  which  seemed 
most  relevant  to  the  immediate  problem. 
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DEFINITION  OF  TERMS 


In  general  the  terminology  of  ultrasonics  and  ordinary 
sonic  phenomena  is  the  same;  since  ultrasonics  is  but  a 
special  branch  of  acoustics*  In  order  to  avoid  confusion, 
the  definition  of  a few  specialized  terms  is  here  included* 
Ultrasonics . The  term  "ultrasonic"  will  be  used  to  denote 
those  frequencies  which  are  beyond  the  limits  of  human  audi- 
bility. These  will  include  all  frequencies  greater  than 
20,000  cycles  per  second  (20  kc).  The  term  "supersonic"  will 
be  reserved  for  phenomena  that  occurs  at  a speed  greater  than 
the  speed  of  sound  in  air.  It  will  be  noticed,  that  in  the 
literature  both  terms  are  used  interchangeably. 

Cavitation.  When  an  ultrasonic  beam  of  high  enough  intensity 
passes  through  a liquid,  bubbles  or  cavities  are  formed  at 
regularly  defined  points  in  the  liquid.  This  phenomena  has 
been  called  cavitation.  In  this  paper  the  term  "cavitation" 
will  be  used  to  describe  both  the  formation  of  gas  bubbles, 
as  well  as  the  formation  of  natural  holes  or  cavities  in  the 
liquid. 

Transducer.  Any  element  which  can  be  mechanically  or  elec- 
trically excited  to  produce  or  receive  ultrasonic  vibrations. 
Group  Velocity.  This  terms  applies  to  the  velocity  with  which 
the  envelope  of  a wave  is  propagated  when  the  wave  is  ampli- 
tude modulated.  This  is  the  most  important  velocity  encountered 


on#  s '.:?tS.Zu  'i  - vr • Ion  ' ~n . d ?r>  r£->  m 

. 

• •••  f-J:  ’•  n . d ’ Jb  9\.  . ,'oi 

. 

. 

•-  : Jfr  dj3  -•  .0.00  , ■ d ^n&  or  ' • 

. 

, 

. £ up IX 

. 

. 

■ t 

n±  asiddvjso 

. 

J ' .. 

- 

- 

. 


iv 


in  ultrasonics.  Unless  otherwise  stated  the  term  velocity 
will  be  used  to  mean  group  velocity. 

Longitudinal  Waves.  The  term  longitudinal  wave  applies  to 
that  wave  motion  where  the  particle  displacement  is  parallel 
to  the  direction  of  propagation  of  the  wave.  Since,  as  a 
rule,  only  longitudinal  waves  can  be  transmitted  through  a 
liquid  medium,  these  will  be  the  most  important  in  this 
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LIST  OF  SYMBOLS  USED 


1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 

9. 


10. 


11. 


12. 


13. 


14. 


~ wavelength 

c - velocity  of  a sound  wave 
y r ratio  of  specific  heats 
S*  s density 
Is  - frequency 

ytx  = velocity  of  vibrating  particle 
p0  a hydrostatic  pressure 
pv  s vapour  pressure 
^ s viscosity 
(p  r surface  tension 
K Boltzmann' s constant 
R = Universal  gas  constant 
T = absolute  temperature 
A s Planck's  constant 
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INTRODUCTION  TO  ULTRASONICS 


(1) 


1.  General  Properties 

Ultrasonic  waves  are  essentially  plane  waves.  Eecause 
of  their  short  wavelength,  rectilinear  propagation  is  char- 
acteristic of  ultrasonic  waves.  Generally,  beams  are  propa- 
gated from  a plane  transducer  with  very  little  divergence.  In 
an  ordinary  quartz  plate  the  energy  is  located  in  a cone  whose 
half  angle  of  spread  is  given  by 


If  the  medium  through  which  the  beam  is  propagated  is 
full  of  small  discontinuities  the  ultrasonic  energy  may  be 
broken  up  and  reflected  and  refracted  haphazardly.  This  is 
an  important  characteristic  in  considering  cavitation  in 
liquids.  In  general  it  is  usual  to  assume  that  the  ultra- 
sonic vibrations  take  place  too  quickly  for  heat  to  exchange. 
Hence  we  may  write  an  expression  for  the  velocity  of  ultra- 
sonic waves  in  liquids  as 


sin  A = 1.2A 


where 


A 5 half  angle  of  spread 


D * distance  across  transducer  face 


C 


(1)  B.  Carlin,  Ultrasonics , McGraw  Hill,  New  York  (1949) 
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where 

/3  is  = isothermal  compressibility 
/3  ad  - adiabatic  compressibility 
The  dissipation  of  energy  in  an  ultrasonic  system  is 
due  primarily  to  the  viscosity  or  damping  of  the  medium. 

There  is  also  some  energy  loss  due  to  heat  conduction  and 
radiation  between  compression  and  rarefaction  sections* 

2.  Range  of  Frequencies  and  Wavelengths 

From  the  familiar  relationship  between  frequency,  , 
wavelength,  /\  and  velocity, c 

c = A > 

it  is  possible  to  calculate  the  frequency  and  wavelength 

ranges  in  gases,  liquids,  and  solids  if  the  velocity  in  that 

( 2) 

medium  is  known.  Eergmann'  ' gives  for  ultrasonic  waves  in 
air  (velocity  of  sound  c = 330  m/s)  a wavelength  range  of  1.6 

cm.  to  0.6  x 10”^cm. , in  liquids  (cc=1200m/s)  a range  of  6 cm. 

-4 

to  2.4  x 10  cm.,  and  in  solids  ( cc?4000  m/s)  a range  from 
~ -4 

20  cm.  to  8 x 10  cm. 

3. Generation  of  Ultrasonic  Waves 

The  method  of  generating  high  frequency  sound  waves 
depends  upon  the  power  output  necessary  and  the  frequency 
range  required.  The  greatest  amount  of  energy  is  commonly 
produced  in  ultrasonics  by  systems  designed  for  agitation 
and  cavitation  experiments.  Energy  outputs  of  10  watts/cm. 

(2)  L.  Bergmann,  Per  Ultraschall,  translated  by  H.  S.  Hat- 
field, Ultrasonics  and  their  Scientific  and  Technical 
Applications,  John  Wiley  and  Sons,  New  York  (1942) 
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have  been  reported. ^ ^ The  most  common  way  of  producing 
ultrasonic  vibrations  is  by  a magnetostriction  rod  or 
piezoelectric  crystal  transducer.  The  former  method  is  used 
when  frequencies  up  to  100  kc.  are  needed  while  the  latter 
method  is  used  for  those  frequencies  in  the  range  from  100  kc. 
to  15  me.  In  recent  years  these  two  methods  have  been  used 
to  the  exclusion  of  others  in  most  cases. 

4.  Detection  of  Ultrasonic  Waves 

When  an  ultrasonic  beam  is  projected  against  a boundary 
it  produces  two  distinct  effects.  First,  there  is  an  alter- 
nating pressure  at  the  frequency  of  wave  propagation,  and 
secondly,  direct  pressure  due  to  the  radiation.  From  any 
text  on  acoustics  it  is  shown  that  the  alternating  pressure 
is  determined  by 

J=  Pm 

3l 

where 

J 3 intensity/cm.2-sec. 

A 3 maximum  amplitude  of 

the  vibrating  particles 
P 3 sound  pressure 
C 3 velocity  of  sound  wave 
of  frequency  w/2tr 

We  see  that  the  measurement  of  sound  intensity  is  possible 
only  if  one  of  the  magnitudes  A, m , or  P can  be  measured  or 
alternatively,  if  both^U  and  P can  be  measured.  In  the  ultra- 
sonic range  it  is  Impossible  to  accomplish  this  with  mechani- 
cal receivers,  since  these  have  dimensions  too  great  compared 


(3)  B.  Carlin,  Ref.  1 
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with  the  sound  waves.  However,  since  we  are  dealing  here  with 
high  frequencies  we  have  also  comparatively  high  intensities 
and  are  enabled  to  make  use  for  measurement,  of  the  second 
effect  mentioned  above,  the  sound  radiation  pressure.  This 
unidirectional  pressure  finds  its  theoretical  explanation  in 
the  nonlinearity  of  the  wave  equation  when  account  is  taken 
of  the  quadratic  terms.  Sound  radiation  pressure,  S,  is  given 
by  the  relationship 

$•  i0r+i)4 

The  measurement  of  sound  radiation  pressure  is  usually  made 
by  radiometers. 

5.  General  Theory  of  Ultrasonic  Wave  Motion 

From  any  text  on  elasticity , we  have  for  the  general 
equation  of  the  displacement  of  a particle  in  any  medium  in  a 
periodic  force  field 

(h  + S)  + & V M - s % ■ ° 

where 


•AL  - displacement  of  the  particle 
in  the  q]_  direction 
A = Lame's  constant 
Gr  = modulus  of  rigidity 
e * volume  expansion 
y - density  of  medium 


This  is  a set  of  differential  equations  of  motion  represent- 
ing types  of  displacement  in  which  the  medium  obeys  Hooke's 


(4)  8.  Timoshenko,  Theory  of  Elasticity,  McG-raw  Hill, 

New  York  (1936T"" 
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law  and  remains  in  a static  state  as  a whole,  i.e.,  the 
only  forces  acting  are  internal  stresses  and  shears  and 
those  due  to  particle  acceleration  forces. 

The  above  represents  a set  of  theoretical  equations  of 
motion  subject  to  physical  and  boundary  conditions.  Since 
longitudinal  waves  are  the  only  type  of  ultrasonic  waves 
possible  in  liquids,  it  is  possible  to  assume  physical  con- 
ditions and  derive  the  general  equation  of  motion  of  the 
particles  of  an  ideal  liquid  for  such  waves.  A longitudinal 
wave  has  been  defined  as  one  in  which  the  particle  velocity 
is  restricted  to  a motion  in  the  direction  of  the  wave  propa 
gation.  Such  waves  are  dilatational  and  irrotational;  i.e., 
there  is  a change  in  the  body  volume  without  rotation.  Now 
since  we  have  the  particle  disolacement M as  the  only  motion 
then  the  volume  expansion  in  the  q^  direction  becomes 


= -2-  is.  = 


assuming  e as  a function  of  alone.  Also,  since  the  motion 
is  restricted  in  the  direction  of  propagation,  if  we  take 
this  as  q-p  then  since 


v 


7 


31 


we  have 


Vm  = 


Hence  our  general  equation  of  motion,  essentially  the  wave 
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equation,  becomes 

(A  -MG--)  = f kkL 

In  actuality,  there  is  always  some  shear  stress  since  the 
edges  of  the  path  cannot  move  back  and  forth  quickly. 

Solutions  of  the  above  differential  equation,  with  the 
application  of  further  physical  and  boundary  conditions,  will 
describe,  theoretically,  an  ultrasonic  wave  field. 

6.  History  of  Ultrasonics 

The  application  of  ultrasonics  in  the  solution  of  in- 
dustrial problems  has  gained  much  popularity  in  the  last 
decade.  Although  phenomena  produced  by  high  frequency  sound 
waves  was  observed  over  half  a century  ago,  the  science  of 
ultrasonics  was  slow  in  its  development.  This  was  due  primar- 
ily to  inadequate  equipment  and  lack  of  knowledge  of  the  basic 
structure  of  matter. 

Prior  to  the  first  world  war  ultrasonic  waves  were  pro- 
duced by  electric  sparks  or  some  mechanical  or  thermal  vibra- 
tor. Langevin,  in  his  work  on  under  water  signalling,  was 
the  first  to  apply  the  piezoelectric  effect  to  the  generation 
of  ultrasonic  waves.  This  provided  the  first  dependable 
source  of  ultrasonic  vibrations  of  controllable  frequency  and 
intensity.  Using  quartz  crystal  he  apolied  the  principle  of  the 
piezoelectric  effect  discovered  by  Pierre  and  Jacques  Curie 
in  1880  to  produce  high  frequency  vibrations  for  the  purpose 
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of  underwater  signalling  and  submarine  detection.  (5) 

Langevin,  during  his  experiments,  noticed  the  lethal  effects 
of  ultrasonic  waves  of  high  intensity  upon  small  marine  life.^ 
This  was  perhaps  the  first  time  that  the  cavitation  effects  of 
ultrasonic  waves  was  mentioned.  Little  was  done  at  the  time 
to  further  investigate  and  explain  its  effects.  The  fact 
that  more  power  was  made  available,  if  the  quartz  transducer 
was  vibrating  at  its  resonant  frequency,  was  also  discovered 
by  Langevin.  In  1922  Cady  boosted  the  advance  of  ultrasonic 
research  when,  as  a result  of  his  studies  of  piezoelectric 
resonators,  he  came  forth  and  developed  the  quartz  frequency 
stabilizer.  Three  years  later,  Pierce  gave  new  impetus  to 
research  by  utilizing  a vacuum  tube  circuit  containing  a 
piezoelectric  plate  for  measurement  of  ultrasonic  velocity  in 
carbon  dioxide. 

In  the  period  1925-1927  Boyle,  and  his  colleagues  in 
Canada,  made  a comprehensive  study  of  the  properties  of  these 
waves;  their  velocity,  reflection,  refraction,  interference 
and  diffraction  in  liquids  was  investigated.  It  was  during 
the  course  of  these  experiments  that  cavitation  in  liquids 
produced  by  an  ultrasonic  beam  was  first  investigated  and 
reported.  ^ ^ (9) 

(5)  J.  C.  Hubbard,  Journ.  Acoust.  Soc.  Am., 4, 99  (1932) 

(6)  B.  Carlin,  Ref.  1 

(7)  R.  W.  Boyle  and  G.  B.  Taylor,  Phys.  Rev.,  (11), 27,518  (1926) 

(8)  G-.  B.  Taylor  and  D.O.  Sproule,  Roy . Soc . Can.  Trans . , (III), 

23,91  (1929) 

(9)  R.  W.  Boyle,  et  al,  Roy. Soc. Can. Trans. ,( III) , 23,187  (1929) 
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The  phenomena  was  observed  for  ultrasonic  waves  where 
the  energy  was  concentrated  in  the  form  of  a narrow  angle 
beam.  The  maximum  sound  energy  which  could  be  transmitted 
through  a liquid  medium  was  found  to  be  limited  by  cavitation. 
If  the  hydrostatic  pressure  p0  is  considered  as  the  only  pre- 
ventive of  disruption  in  the  medium,  then  the  maximum  energy 
transmi table  per  square  centimeter  is  given  by 


Experiments  showed  that  cavitation  is  more  readily  produced 
in  highly  volatile  liquids,  hence  it  was  suggested  that  the 
vapour  pressure,  Py,  be  taken  into  consideration.  This  gives 
for  the  maximum  energy  transmi table 

/cm*  = (-ft- 

Experiments  showed  that  the  cavitation  point  was  far  below 

that  given  by  the  above  equation.  It  is  now  known  that  other 

factors  besides  those  listed  above  limit  the  power  output. 

Study  of  the  effects  of  intense  ultrasonic  radiation 

(10) 

began  in  1927  with  Wood  and  Loomis  who  carried  out  exten- 
sive experiments  using  a specially  designed  vacuum  tube  gener- 
ator capable  of  2000  watt  output. Shortly  thereafter  it 
became  apparent  that  emulsification,  the  destruction  of 
bacteria,  the  freeing  of  antibodies  from  bacterial  cells, 


(10)  R.W.Wood  and  A. L. Loomis,  Phys.  Rev.,  11,  29  (1927) 

(11)  R.W.Wood  and  A. L. Loomis,  Phil. Mag.,  4,  417  (1927) 
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the  acceleration  of  chemical  reactions,  the  production  of 
luminescence  and  several  other  phenomena,  observed  in  an 
intense  ultrasonic  wave  field,  were  dependent  for  their 
existence  upon  cavitation.  It  was  also  observed  that  in  de- 
gassed liquids,  or  in  liquids  under  high  hydrostatic  pressures, 
most  of  these  effects  either  disappeared  or  proceeded  very 
feebly. 

( 12) 

In  1932  Debye  and  Sears  performed  light  scattering 
experiments  with  ultrasonic  waves  in  which  they  observed 
diffraction  patterns.  There  is  a possibility  that  this  pheno- 
mena is  due  to  cavitation.  In  the  early  1930's,  physicists 
and  chemists  used  ultrasonics  to  further  study  specific  heats 
of  gases  with  respect  to  the  relative  parts  played  by  vibrat- 
ional and  rotational  frequencies  of  molecules  at  high  tempera- 
ture. This  was  approached  through  a study  of  absorption  and 
velocity  determination  in  gases  as  well  as  in  liquids.  Much 
research  along  these  lines  is  still  in  progress. 

Further  progress  in  the  investigation  of  ultrasonic 
cavitation  had  to  wait  for  a new  theory  of  the  physical 
structure  of  liquids.  H.  Eyring  in  1936  advanced  some  new 
ideas  on  liquid  structure  which  was  extensively  developed  by 
R.  Furth  in  1941.  This  theory  will  be  discussed  more  thor- 
oughly later. 

The  greatest  avenues  of  further  development  of  ultra- 
sonic research  lie  in  the  fields  of  metallurgy  and  biophysics. 

(12)  P.  Debye  and  F.W. Sears,  Proc.Natl.  Acad.  Sci.U.S., 

18,410  (1932) 
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Before  any  real  progress  can  evolve  in  these  fields,  a 
concise  mathematical  and  physical  analysis  of  ultrasonic 
cavitation  is  necessary.  A great  deal  remains  to  be  accom- 
plished 
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EXPERIMENTAL  PROCEDURES  IN  THE  INVESTIGATION 
OF  ULTRASONIC  CAVITATION^ 13 ^ 

1.  Apparatus 

The  following  is  a brief  account  of  the  experimental 
procedure  and  apparatus  used  by  Brlggs^-^^  and  his  colleagues 
during  a series  of  experiments  designed  to  test  the  cavita- 
tion strength  of  various  liquids.  A sectional  view  of  the 
test  chamber  and  crystal  units  used  for  studying  cavitation 
is  shown  in  Fig.  it  The  major  faces  of  the  crystals  were 
coated  with  gold  and  cemented  together  in  groups  of  four. 

This  crystal  group  was  then  cemented  to  a quarter  wave 
resonator  block  which  made  a half  wave  unit  with  the  node  at 
the  interface  between  the  crystal  and  resonator  block.  This 
unit  fitted  into  a cup  faced  with  rho-c  rubber  (a  rubber 
having  the  same  density  and  dilation  velocity  as  water), into 
which  was  placed  the  liquid  to  be  tested* 

The  radiator  radiated  waves  through  the  liquid,  rubber 
and  standpipe  to  a second  crystal  unit  where  they  were  picked 
up  and  generated  into  an  electrical  replica  which  was  put  on 
a cathode  ray  tube.  In  order  to  avoid  cavitation  in  the 
standpipe,  it  was  filled  with  degassed  castor  oil.  By  the 
use  of  a vacuum  pump  and  a rotary  oil  pump,  the  liquid  to  be 

(13)  H. B. Briggs, et  al,  Journ.  Acoust .Soc.Am. , 19,664  (1947) 

* All  figures  and  tables  mentioned  may  be  found  in  Ref.  13 
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tested  for  cavitation  resistance  was  also  thoroughly 
degassed. 

A schematic  diagram  of  the  power  source  and  circuits 
used  for  testing  liquids  is  shown  in  Fig.  2.  The  power 
input,  W,  to  the  quartz  crystal  projector  was  calculated  from 
the  voltage  V and  resonance  resistance  R of  the  projector 
using  the  relationship 

W = V2/R 

Two  power  supplies  were  available;  one  a 75  watt  supply  used 
for  either  pulsing  or  steady  power,  the  other  a 1000  watt 
supply  used  only  for  pulsing  power. 

For  steady  power  input,  the  voltage  across  the  crystal 
terminals  was  read  off  directly  on  a vacuum  tube  voltmeter, 
while  the  minimum  resistance  was  taken  as  the  resonance 
resistance.  For  pulsing  tests  the  voltage  applied  to  the 
projector  was  determined  from  the  amplitude  of  the  oscillo- 
scope deflection.  The  oscilloscope  also  served  to  measure 
the  relative  magnitude  of  response  picked  up  by  the  receiver 
for  varying  inputs.  It  was,  therefore,  possible  to  draw 
input-output  curves  as  shown  in  Fig.  3 for  non-degassed  kero- 
sene. Wave  form  or  pulse  form,  as  shown  by  the  oscilloscope, 
indicated  when  cavitation  was  present. 

2.  Procedure 

Assuming  the  radiator  to  be  100  percent  efficient, 
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peak  pressures 
formula 


or 


can  be  calculated  approximately  by  the 

w - 'a  fJ  A/st 

k.  ■ 


where 

pm  * peak  pressure  in  dynes/cm2 
A - radiating  area  in  cm2 
W - power  output  in  ergs/sec 

As  an  example,  we  examine  the  input-output  curve  of 

Fig.  3*  The  relationship  between  input  and  output  is  linear 

up  to  800  volts  (cavitation  point)  after  which  it  bends  and 

rises  very  slowly.  It  was  found  that  the  resonance  resist- 
ed 

ance  was  2x10"^  ohms.  With  a projector  surface  area  of 

p 

6.45  cm  we  get  for  the  power  output  per  unit  area  of  projector 


W/cm2  * V2/RA  = 0.5  watts/cm2 


If  the  density  of  non-degassed  kerosene  is  taken  as  0.81 
gm/cm2  and  the  velocity  of  the  wave  in  the  liquid  as  1.324  x 
10^  cm/sec,  we  get  for  the  maximum  pressure  a value  a 
little  larger  than  10^  dynes/cm2,  which  is  roughly  atmos- 
pheric pressure.  Hence,  it  would  seem  safe  to  conclude  in 
this  case  that  a light  undegassed  liquid  is  held  together 
principally  by  atmospheric  pressure  and  does  not  exhibit  any 
cohesive  force. 


3.  Results  of  the  Experiments 

For  a better  understanding  of  the  process  of  cavitation, 
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Briggs  made  measurements  on  liquids  having  wide  ranges 
of  viscosity,  vapor  pressure  and  chemical  structure.  Table 
I lists  some  of  the  properties  of  these  liquids.  It  is  to 
be  noted  that  the  vapour  pressures  are  not  given  specifi- 
cally but  rather  as  lying  within  certain  ranges. 

The  first  series  of  tests  were  made  on  the  effects  of 
dissolved  gases  on  cavitation  voltage.  By  varying  the 
pressure  on  the  liquid  the  effect  of  hydrostatic  pressure  was 
determined.  Fig.  4 shows  the  voltage  across  the  crystal  at 
which  cavitation  occurred  as  a function  of  the  ambient 
pressure  for  degassed  kerosene.  The  resultant  curve  is  a 
straight  line  which  when  extended  cuts  the  axis  at  a negative 
pressure  of  60  cm.  of  mercury.  Hence  degassed  kerosene  will 
stand  a negative  pressure  of  nearly  one  atmosphere  before 
cavitation  occurs.  We  can  also  conclude  that  the  cavitation 
pressure  is  the  sum  of  the  cohesive  pressure  (tensile  strength) 
plus  the  hydrostatic  pressure. 

Measurements  which  were  made  on  the  cohesive  pressure 
of  the  liquids  used  are  given  in  the  last  column  of  Table  I. 

Pc  is  the  peak  acoustic  pressure  at  which  cavitation  begins 
at  one  atmosphere  hydrostatic  pressure  in  the  liquid.  The 
cohesive  pressure  p_  is  then 

Pc  ' PC  - 1 

The  value  of  p depends  to  a large  extent  on  the  previous 
c 

(14)  H.  B.  Briggs,  et  al,  Ref.  13 
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history  of  the  liquid.  A cavitated  liquid  is  "weakened" 
and  will  require  some  time  to  return  to  its  earlier  state; 
if  the  experiment  is  repeated  the  liquid  will  cavitate  at  a 
lower  acoustic  pressure* 

It  was  found  that  the  only  quantity  consistently  re- 
lated to  the  cohesive  pressure  is  the  viscosity.  Fig.iT 
shows  cohesive  pressure  as  a function  of  viscosity,  the  latter 
being  plotted  on  a logarithmic  scale.  Briggs  found  the 
equation  of  this  curve,  which  is  a straight  line  within  ex- 
perimental error,  to  be 

^ ( Vt,)  • 2.0S-& 

where 

h 

r = viscosity  in  poises 
• © * 0.0013  poise 

pc  i cohesive  pressure  in  atmospheres 

From  the  above  equation  it  is  seen  that  as  h.  approaches 
*l0  the  cohesive  pressure  approaches  zero.  Hence  when  the  vis- 
cosity approaches  0.0013  poise  the  liquid  is  held  together 
only  by  the  static  pressure  of  the  atmosphere.  The  theoreti- 
cal interpretation  of  the  cavitation  point  of  a liquid  as  a 
function  of  viscosity  will  be  given  when  Eyring's  theory  of 
viscosity  is  considered. 
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PRODUCTION  AND  PROPERTIES  OF  CAVITIES  IN  LIQUIDS 
1.  Production 

When  a piezoelectric  quartz  crystal  is  immersed  in  a 
liquid  and  made  to  vibrate  at  ultrasonic  frequencies,  bubble® 
or  cavities  are  observed  in  the  liquid  if  conditions  are  fav- 
orable. The  simplest  arrangement  is  probably  found  in  the 

ultrasonic  interferometer,  a description  of  which  is  found  in 
( 15) 

Colby.  If  a solid  reflecting  surface  is  placed  an  integral 

number  of  high  wave  lengths  from  the  transducer,  a stationary 
wave  is  set  up  in  the  liquid  with  a displacement  node  at  the 
face  of  the  transducer  and  also  at  the  reflecting  wall.  Since 
in  a stationary  wave  the  displacement  wave  and  compression  wave 
are  ninety  degrees  out  of  phase,  we  find  pressure  loops  and 
displacement  nodes  coinciding.  Cavities  are  formed  during 
the  negative  half  of  the  pressure  cycle,  and  are  observed  to 
collect  at  pressure  loops  of  the  stationary  wave  pattern. 

Hence,  we  find  cavities  present  at  the  faces  of  the  transducer 
and  reflector  and  also  at  distances  of  ^/2,  where  A is  the 
wavelength,  from  the  transducer.  It  will  be  sufficient  at 
this  point  to  explain  the  formation  of  these  cavities  as  a 
"tearing  apart"  of  the  molecules  of  the  liquid,  which  leaves 
empty  spaces  behind. 

(15)  M.  Y.  Colby,  Sound  Waves  and  Acoustics,  Henry  Holt  and 
Co.,  New  York  Tl93S) 
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2.  Properties 

According  to  Sollner^1^^  these  cavities  may  be  made 
visible  by  irradiating  long  layers  of  completely  degassed 
liquid  with  an  ultrasonic  beam  of  high  intensity.  Due  to  the 
nearly  sinusoidal  pressure  distribution  in  the  liquid,  we  have 
zones  of  greater  density,  therefore  less  transparency  to 
light.  These  zones  disappear,  when  the  beam  is  shut  off, 
without  leaving  bubbles  behind.  An  extension  of  this  phen- 
omena is  used  in  the  design  of  the  ultrasonic  diffraction 
"grating". 

An  ingenious  method  of  determining  experimentally  the 
loci  of  equal  phase  in  an  ultrasonic  compressional  wave  beam 
without  utilizing  stationary  wave  patterns,  has  been  devised 
by  Labow.^^)  The  electric  signal  generated  by  the  piezo- 
electric action  of  a Rochelle  salt  microphone,  which  has  the 
same  phase  as  the  compression  wave  in  the  liquid,  is  combined 
with  a small  electric  signal  of  the  same  frequency  having 
the  phase  of  the  transducer  voltage.  The  result  of  this 
superposition  is  a sinusoidal  spatial  variation  of  the 
signal  generated  by  the  microphone  as  the  latter  is  moved 
in  the  direction  of  the  beam.  The  wave  front  is  traced  by 
following  a maximum  of  this  sinusoidal  variation  as  the 

(16)  K.  Sollner,  Trans.  Faraday  Soc.  32,  1537  (1936) 

(17)  L.W.Labow,  Journ.  Acoust.  Soc.  Am.,  17,19  (1945) 
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microphone  is  made  to  move  perpendicular  to  the  direction  of 
propagation  of  the  beam.  Hence  it  is  possible  to  determine 
where  cavitation  should  take  olace  if  it  is  not  visible. 
However,  with  proper  power  input,  frequency  and  a suffic- 
iently gassed  liquid,  the  vapor  dissolved  in  the  liquid  soon 
fills  up  the  holes  caused  by  the  tearing  apart  of  the  mole- 
cules, forming  gas  bubbles  which  grow  until  sufficiently 
big  to  be  visible. 

These  spherical  bubbles  of  gas  pulsate  radially  under 
the  influence  of  the  alternating  sound  pressure. The 
pulsating  bubble  has  potential  energy  from  the  elasticity  of 
the  gas  and  kinetic  energy  arising  from  the  movement  of  the 
liquid  in  the  vicinity.  We  have  resonance  occurring  at  that 
frequency  when  the  mean  potential  and  kinetic  energies  are 
equal.  Smith  has  shown  that  the  lowest  resonance  frequency 
of  a bubble  formed  by  cavitation  is  given  by 


For  water  under  standard  conditions  a bubble  of  0.033  mm* 
radius  has  a resonant  frequency  of  100  kc.  These  resonant 


where 


f - resonant  frequency 
r - radius  of  the  bubble 

(T  - surface  tension  of  liquid  gas  interface 


(18)  F.D. Smith,  Phil. Mag.,  19,  1147  (1935) 
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bubbles  are  characterized  by  violent  activity  and  appear  to 
the  eye  as  a luminous  streak  on  account  of  the  persistence 
of  vision.  Larger  and  smaller  bubbles  are  static,  since  due 
to  a change  in  radius  the  bubbles  no  longer  resonate  to  the 
frequency  of  the  applied  field.  "The  irregular  motion  of  the 
bubbles  in  such  a field  may  possibly  be  explained  as  due  to 
Bjerknes  forces  set  up  between  the  bubbles  and  the  sound 
field  with  attraction  or  repulsion  between  the  bubbles  de- 
pending on  the  relative  phase  of  their  vibrations.  The  motions 
of  the  bubbles  bear  some  similarity  to  the  behavior  of 
charged  particles  in  a varying  electrostatic  field.  "^9) 
Richards ^ 20 ^ has  observed  that  cavitation  is  accompanied  by 
severe  electrical  disturbances  in  liquids.  Over  distances 
molecular  in  size  these  may  be  great  enough  to  produce  ioni- 
zation, combination,  etc*  which  may  explain  the  behavior  of 
these  gas  bubbles. 

Some  of  the  properties  of  these  bubbles,  formed  by 
cavitation,  are  brought  out  if  we  follow  the  derivation  for 
the  frequency  of  a resonant  bubble  given  above  by  Smith. 

Bubbles  are  small  in  comparison  with  the  wavelength,  hence  the 
alternating  sound  pressure  acts  uniformly  over  the  surface 
of  the  bubble  and  they  pulsate  in  the  simplest  radial  mode. 

Each  bubble  has  a restoring  force,  inertia  and  damping  and 
pulsates  as  a resonant  mechanical  system  which  has  one  degree 
of  freedom. 

(19)  H.B. Briggs,  etal,  Joum.Acoust.  Soc .Am. , 19,664  (1947) 

(20)  W.T. Richards,  Journ.Appl.Phys. , 9,  298  (1938) 
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The  restoring  force  per  unit  area  which  arises  from  the 
elasticity  of  the  gas  in  the  bubble  is 

r - X(r.i-v) 

We  will  assume  here  that  the  compressions  and  expansions 
of  the  gas  in  the  bubbles  take  place  too  fast  for  any  heat 
transfer,  hence  an  adiabatic  process.  In  the  last  expression 
Pis  the  restoring  force  per  unit  area,  while  ,r,  p0  and  (T 
have  the  same  meaning  as  before. 

The  effective  mass  per  unit  area,  oC  , of  a pulsating 
bubble  arises  from  the  kinetic  energy  of  the  liquid  moved 
radially  by  the  pulsations^ and  is  given  by 

= ? A.  ( I t /?aA3) 

whe  re 


The  radiation  resistance  or  damping  factor  per  unit 
area  /3  arises  from  the  scattering  of  the  sound  by  the 
bubble^22)  and  is  given  by 

a = f c a*v/(,  + *v) 

where 

c - velocity  of  sound  in  the  liquid 

(21)  R.W.Wood  and  A. L. Loomis,  Phil. Mag.,  4,417  (1927) 

(22)  J.B. Crandall,  The  Theory  of  Vibrating  Systems  and  Sound 

D.  Van  Nostrand  and  Co.,  New  York  (1926) 
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Now, if  we  call  p equal  to  the  instantaneous  alternating 
pressure  in  the  liquid  in  the  absence  of  a bubble,  then  the 
equation  of  motion  of  the  bubble  is  given  by 


dti 

cLt* 


where  § is  the  displacement  of  the  surface  of  the  bubble 
from  its  equilibrium  position.  This  is  analogous  to  elec- 
trical and  mechanical  forced  oscillatory  systems. 

From  the  theory  of  forced  oscillations  we  have  resonance 
occurring  when 


Ma 


+ rt 


0 


or  when 


That  is, 
given  by 


resonance  will  occur  when  we  have  a frequency 


I 

c37 T 


'/a 


f 


If  we  substitute  for  f**  and  we  have 

^ ^ — r 3 ir ( pa  + 3. (r/\}  ~i 

T l $(  /V  J 


/a 


which  gives  the  resonant  frequency.  In  his  derivation  for 
^ Smith  assumed  that  for  water  k2r^<(l,  hence  his  formula 


follows.  It  may  well  be  that  these  resonant  bubbles  are  the 
cause  of  the  destructive  phenomena  associated  with  cavitation 
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STRUCTURE  OF  LIQUIDS  AND  THEIR  PROPERTIES 
IN  AN  ULTRASONIC  FIELD 


Furthv  J has  suggested  a new  theory  of  liquid  structure 
in  which  holes  in  a continuum  are  treated  by  the  methods  of 
statistical  mechanics . ^ This  is  exactly  the  reverse  of 
the  treatment  of  solid  particles  in  empty  space  encountered 
in  kinetic  theory.  These  holes  are  assumed  to  contain 
vapour  at  the  vapour  pressure  corresponding  to  the  tempera- 
ture. This  differs  from  Eyring's^2^)  theory  of  holes  in 
that  Earing  assumes  the  existence  of  natural  holes  in  a 
liquid,  which,  when  in  the  equilibrium  state,  are  devoid  of 
vapour. 

Furth's  treatment  assumes  the  existence  of  a known 
metastable  state  in  which  the  temperature  exceeds  the  satura- 
tion value  corresponding  to  the  pressure;  the  vapour  pressure 
exceeding  the  hydrostatic  pressure.  For  the  limit  of  this 
metastable  state  he  obtains  the  following  equation 

_ jp  LI  c r3/a_ 
r/yr  /o  ~ (/<r)//at 

where 

T - temperature 


(23)  R.Furth,  Proc.  Camb.  Phil.  Soc.,  37,252  (1941) 

(24)  R. 8. Silver,  Nature,  150,605  (1942) 

(25)  H.Eyring,  J.Chem.,  Phys.  4,283  (1936) 
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For  example,  for  water  at  normal  temperature  we  get 
Pv  - pQ  * 8x10^  lb s/in2. 

Now  if  pQ  is  considered  atmospheric  pressure  then  the 
excess  of  vapour  pressure  demands  a high  temperature  before 
it  can  be  realized,  and  in  this  case  has  little  meaning 
since  the  pressure  is  far  above  the  critical  breaking 
pressure  of  the  liquid.  In  reality  the  right  hand  side  of 
the  above  equation  needs  to  be  altered  to  allow  for  differ- 
ent values  of  the  temperature.  As  the  equation  stands  above 
the  two  sides  become  equated  for  water  at  a temperature  of 
about  337°C  when  the  vapour  pressure  is  in  excess  of 
2000  lbs/in2. 

The  superheating  of  water  has  never  been  attained  to 
this  degree,  due  to  the  presence  of  particles  and  dissolved 
gasses.  The  highest  reached  is  about  180°C  for  water  under 
artificial  conditions. 

There  is  another  way  of  interpreting  this  equation.  It 
is  usually  assumed  that  the  hydrostatic  pressure  is  equal  to 
or  less  than  zero,  whence  at  a given  temperature  pv  - pQ  is 
less  than  or  equal  to  the  vapour  pressure.  However,  if  the 
hydrostatic  pressure  is  taken  less  than  zero  then  the  relation- 
ship between  p^  - pQ  and  the  vapour  pressure  given  above  is  no 
longer  true.  From  the  above  data  it  appears  that  the  equation 
defining  the  limit  of  the  metastable  state  would  be  satisfied 
if  the  hydrostatic  pressure  were  of  the  order  of  8 x 10^1bs/in2, 
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since  the  vapour  pressure  at  standard  temperature  is  low 
enough  to  be  negligible.  From  this  we  might  infer  that  the 
liquid  should  theoretically  have  a tensile  strength  up  to 
the  limit  of  the  hydrostatic  pressure.  However,  particles 
and  dissolved  gases  in  the  liquid  tend  to  lower  the  limit 
far  below  the  expected  value.  From  elementary  considerations 
it  can  be  shown  that  in  the  case  of  water,  for  example,  the 
tensile  strength  is  negligible.  This  seems  to  be  borne  out 
by  experimental  facts  as  was  shown  earlier  by  Briggs  for  un- 
degassed kerosene.  If  we  employ  the  equation  for  equili- 
brium of  a bubble  of  radius  r,  it  is  found  that  the  tensile 
strength  limit  lies  between  2 <T/r  and  4 O'/Jr,  Bubbles  with 
radii  as  small  as  lCT^cm.  would  make  this  limit  of  the  order 
of  about  one  atmosphere.  Under  ordinary  circumstances  evapo- 
ration nuclei  of  much  greater  radius  exist  in  water,  de- 
creasing still  further  the  tensile  strength. 

/ of) 

Silver'  suggests  the  following  interpretation  to  the 
seeming  fallacy  in  Furth’s  equation  for  the  limit  of  a meta- 
stable state  of  a liquid  in  which  vapour  filled  holes  are 
assumed.  While  evaporation  nuclei  may  exist  in  a liquid, 
evaporation  will  not  take  place  unless  there  exists  a suit- 
able temperature  distribution  formed  to  allow  the  conduction 
of  heat  into  these  nuclei.  Also,  this  must  occur  at  a 
limited  rate  governed  by  the  distribution  of  the  temperature 
and  the  surface  area  of  these  nuclei.  Hence,  we  may  deduce 


(26)  R. 5. Silver,  Nature  150  , 605  (1942) 
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that  the  tensile  stress  above  the  limit  imposed  by  the 
nuclei  may  persist  for  a short  time.  The  value  attained 
depends  on  the  relative  times  for  evaporation  and  also  for 
the  growth  of  the  stress.  The  suggestion  is  that  the  ulti- 
mate limit  corresponding  to  the  hydrostatic  pressure  derived 
from  Furth  could  be  carried  for  an  infinitesimal  time  inter- 
val and  still  be  transmissable  as  a tensile  wave  of  rect- 
angular wave  front.  For  example,  a very  steep  tensile  wave 
front  is  to  be  expected  for  water. 

The  problem  then  is  to  find  a peak  tension  before  the 
liquid  ruptures^ . Consider  a liquid  saturated  with  gas. 

If  some  of  this  liquid  is  enclosed  in  a container  and  tension 
is  applied,  we  have  supersaturation  due  to  the  reduced  pressure. 
A bubble  of  gas  would  form  and  a new  equilibrium  state  would 
be  reached.  Yet  some  time  is  necessary  for  diffusion  to  the 
weakest  boundary  to  form  the  initial  bubble  and  the  saturated 
liquid  would  stand  a momentary  tension.  However,  if  the 
liquid  is  completely  degassed  the  pressure  can  still  be  re- 
duced below  its  own  vapour  pressure  without  rupture,  since 
completely  degassed  liquids  have  a tensile  strength  which 
permits  sound  pressure  higher  than  the  hydrostatic  pressure 
to  be  transmitted. 

If  we  consider  Van  der  Waals  equation  of  state  which 


(27)  H.B. Briggs,  et  al.,  Ref.  13 
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holds  approximately  for  liquids 

(p  + a/v2)(v  - b)  - RT 

we  notice  that  as  the  pressure  p is  varied  the  volume  v 
changes  but  little,  hence  the  term  a/v2  may  be  considered  a 
constant.  It  is  in  the  form  of  an  external  force  that  holds 
the  volume  of  the  liquid  constant  although  its  shape  may 
change.  The  magnitude  of  the  internal  tension  may  be  esti- 
mated from  the  compressibility  of  liquids  which  is  of  the 
order  of  about  Av/v  * 10**^  per  atmosphere  pressure.  The 
Inward  pressure  a/v2  is  of  the  order  of  10^  atmospheres, 
and  application  of  ordinary  pressures  externally  changes  the 
volume  of  the  liquid  but  little.  Attempts  to  measure  the 
static  breaking  strength  of  liquids  give  erratic  values  far 
below  the  theoretical  values. 

Vincent^ made  a study  of  the  stability  of  a gas 
bubble  in  a liquid.  Any  gas  bubble  should  be  unstable  since 
it  is  buoyant  and  will  tend  to  rise  and  enlarge  in  a liquid. 

A temporary  stability  may  be  established  if  the  pressure  inside 
the  bubble  is  greater  than  the  hydrostatic  pressure.  It  is 
clear  that  the  smaller  the  bubble  the  greater  the  internal  * 
pressure  must  be  for  stability.  It  is  fairly  simple  to  cal- 
culate what  the  radius  of  the  largest  bubble  can  be  if  we  let 

(28)  R.S. Vincent,  Pro c. Phy s. Soc. London  53,126  (1941) 

(29)  R.S. Vincent,  Proc .Phy s. Soc .London  55,41  (1943) 
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and 


r » radius  of  the  bubble 
p - internal  pressure. 


Clearly  the  force  of  pressure  on  an  equatorial  cross  section 
of  bubble  is 

Fe  - (p  - P0)"r2 

and  this  is  balanced  (for  stability)  by  the  force  of  surface 
tension  along  the  periphery 

Fp  3 27TrO“ 

The  total  outward  force  is  then 


or 


F - Fe  - Fp 

F - (p  - p0)*r2  - 2"rcr 


Now  if  the  bubble  is  to  be  stable  the  total  outward 
force  must  equal  zero  and  hence 

2<r 

r ~ (p  - P0) 


For  instance,  Vincent  found  that  the  highest  breaking 
strength  in  oil  is  about  7.8  atmospheres.  We  can  now  cal- 
culate the  radius  of  the  largest  stable  bubble  present.  The 
surface  tension  for  oil  is  about  30  dynes/cm.  and  the  hydro- 
static  pressure  is  -7.8  x 10  dynes/cm  . Assuming  that  the 
hydrostatic  pressure  in  the  bubble  is  negligible,  the  bubble 

which  became  unstable  and  caused  the  break  was  no  greater 

—A 

than  7*7  x 10  cm.  just  before  the  break. 

A bubble  of  such  dimensions  has  very  little  buoyancy 
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and  is  eliminated  from  the  liquid  with  difficulty  even  if 
near  the  surface.  If  as  has  been  assumed  by  Furth  and 
Eyring,  mobile  holes  exist  in  a liquid,  then  the  process  of 
rupture  consists  in  one  of  the  larger  holes  becoming  un- 
stable under  the  tension  applied  to  the  liquid  and 
the  gas  that  apoears  arrives  in  the  enlarged  cavity  by 
subsequent  diffusion  from  the  liquid. 

It  can  be  concluded  then  that  the  formation  of  cavities 
of  unstable  size  in  a degassed  liquid  has  to  proceed  unaided 
by  gas  pressure.  At  a very  early  stage  of  development  of 
the  cavity,  externally  applied  tension  on  the  liquid  fails 
to  be  supported  by  any  internal  tension;  the  cavity  then 
fails  to  be  an  abnormal  point  in  the  liquid  and  becomes  a 
discrete  void  enclosed  by  a boundary  where  the  properties 
of  surface  tension  appear.  Now  if  the  external  tension  is 
released  the  cavity  should  collapse.  However,  if  conditions 
are  not  favorable  for  collapse,  a time  lag  is  involved 
before  rupture  occurs.  This  is  a normal  condition  in  all 
the  dynamic  measurements  involved  in  ultrasonic  cavitation 
experiments  made  by  Briggs and  his  colleagues. 


(30)  H.B. Briggs,  et  al.  Ref.  13 
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HYDRODYNAMIC  THEORY  OF  CAVITATION 
1.  Introduction 

A hydrodynamic  theory  of  cavitation  in  liquids  has  been 

( "51 ) 

suggested  by  Kennard  ' for  idealized  liquids.  It  is  well 
to  keep  in  mind  at  this  point  that  the  primary  effect  of  the 
passage  of  an  ultrasonic  beam  through  a liquid  is  a temperoral 
and  spatial  redistribution  of  the  pressure  within  the  liquid. 

The  liquid  will  cavitate  when  its  pressure  sinks  to  a break- 
ing pressure  p^,  then  the  pressure  rises  at  once  to  a fixed 
cavity  pressure  pc  where  Pb£:pc.  Hence  p^  is  the  minimum 
pressure  that  a liquid  can  stand.  Certain  assumptions  need 
to  be  made  here  in  order  to  keep  the  nature  of  cavitation 
from  becoming  too  complex  to  handle  mathematically.  "I  have 
been  told  by  every  mathematical  physicist  I know  that  the 
analysis  of  cavitation  is  a task  beyond  the  ability  of  present 
day  mathematics."^2) 

It  will  be  assumed  here  that  the  pressure  in  the  bubbles 
is  constant,  the  bubbles  are  small  and  contain  only  vapour 
of  the  liquid  when  the  cavity  pressure  equals  the  vapour 
pressure,  and  that  the  motion  is  irrotational  in  a homo- 
geneous liquid  of  nearly  constant  density.  Hence,  the  variation 
of  density  due  to  the  passage  of  an  ultrasonic  beam  can  be 
neglected  here. 


(31)  E.H. Kennard,  Phys.Rev.  63,  172  (1943) 

(32)  W.T. Richards,  Ref. 20 
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2.  Plane  Breaking  Fronts 


The  one  dimensional  equations  of  motion  for  acoustic 
phenomena  are  given  as 

- _.  / q)  >0  ' ^ P _ __  ^ 

o>t  " s > T?  ” y 

where 

p - excess  of  pressure  above  hydrostatic 
t - time 

S * undisturbed  density 

The  derivatives  in  the  above  equations  are  assumed  to  be 
continuous  functions.  These  equations  may  be  applied  gener- 
ally to  the  relative  motion  of  the  liquid  in  the  neighbor- 
hood of  any  point  if  uniformly  moving  axes  are  employed 

relative  to  which  the  velocity  of  the  liquid  near  the  point 

is  very  small. 

It  has  been  already  asserted  that  cavitation  begins 
when  the  pressure  sinks  to  the  breaking  pressure,  hence  when 


p 5 Pfc  then 


< ° 


> 0 


From  the  point  where  p equals  p-D  a definite  breaking  front 
will  advance  toward  both  sides;  that  front  advancing  toward 
the  right  advances  into  regions  where  ^ o while  the 

front  advancing  to  the  left  advances  into  regions  where 
± 0.  Just  ahead  of  each  breaking  front  conditions 
that  held  at  p equals  p^  must  also  hold. 
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In  order  to  find  the  rate  of  advance  of  a breaking 
front  we  notice  what  happens  when  the  front  moves  to  the 
right  a distance  dx  in  a time  dt. 
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From  the  above  diagram  we  see  that  if  is  the  velo- 
city of  the  plane  breaking  front,  then  at  the  beginning  of 
dt  the  pressure  in  the  uncavitated  liquid  is  p^  over  the 
left  hand  face  of  the  layer  F-j_,  while  it  is  p^  *(!£)  dx 
over  the  right  hand  face  of  F2.  At  the  end  of  dt  the 
pressure  over  the  right  hand  face  of  F2  is  also  p^.  Then 
the  change  in  pressure  at  F2  during  the  time  interval  dt  is 
-^jf!^dx,  or  since  is  the  speed  of  the  plane  breaking  front 
relative  to  the  liquid  so  that  dx  equals  U-^dt,  then  the  change 
may  be  expressed  as 


then 
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We  get  then  for  the  velocity  of  the  front 


From  the  diagram  it  is  seen  that  at  the  boundary  of  the 
front  there  is  a discontinuity  in  pressure.  To  find  the 
impulsive  effect  consider  the  motion  of  a layer  of  liquid  dx. 
On  the  right  hand  face  F2  the  pressure  is  nearly  equal  to  the 
breaking  pressure.  It  sinks  to  this  value  as  the  front  comes 
up  to  the  face,  whereas  on  the  left  hand  face  we  have  acting 
the  steady  cavity  pressure.  The  momentum  change  in  the  layer 
dx  during  the  time  dt  is  then  given  by  (pc  - p^Jdt  per  unit 
area.  If  is  equal  to  the  particle  velocity  just  ahead 

of  the  front  and  M,c  is  the  particle  velocity  just  behind  it, 
then  the  velocity  changes  Impulsively  from  to  M <•  as 

the  front  passes,  and  since  we  have  dx  s U^dt  then 


This  gives  us  the  particle  velocity  just  behind  the  plane 
breaking  front.  'The  same  can  be  shown  to  be  true  for  the 
breaking  front  advancing  to  the  left.  It  can  also  be  shown 
that  the  motion  of  the  liquid  as  given  by above  is  hydro- 
dynamically  possible.  There  are  two  basic  arguments  for  the 
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proof  of  this  - one  where  pc  s Pb>  the  other  where 
Pc  y Pb * Both  of  these  arguments  lead  to  the  same  result; 
namely,  that  the  square  of  the  velocity  of  the  breaking 
front  is  equal  to  or  greater  than  the  square  of  the  velocity 
of  the  sound  wave  in  the  liquid.  From  this  it  follows  that 
Ub  - *c,  but  if  this  is  true  for  any  finite  time  interval 
then  cavitation  does  not  really  occur.  We  can  conclude,  then, 
that  a true  breaking  front  always  travels  faster  than  sound 
in  the  liquid.  The  conditions  which  must  exist  for  the  pro- 
pagation of  a plane  breaking  front  may  be  summarized  as 
follows : 


for  propagation  to  the  left. 

All  of  these  derivatives  are  evaluated  just  ahead  of 
the  front  since  no  hydrodynamic  influence  of  the  first  order 
can  be  propagated  through  a liquid  at  a speed  faster  than 
sound.  Once  started,  a breaking  front  will  continue  advanc- 
ing until  it  arrives  at  a point  beyond  which  the  necessary 
conditions  stated  above  are  no  longer  satisfied. 

3.  The  Cavitated  Region 

In  the  cavitated  region  the  pressure  is  uniformly  pc, 


and  also 


^ ft  for  propagation  to  the  right 


= 0 and  the  particle  velocity  , as 


hence  we  have 
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derived  above,  remains  as  it  was  when  the  breaking  front 
passed  by.  Now  let  <p  equal  the  fraction  of  the  space 
occupied  by  bubbles.  Then  if  pQ  = pb,  p : 0 at  the  trail- 
ing edge  of  the  front.  But  if  pc  > p^  then  p starts  from  <po 
which  is  equal  to  the  difference  between  the  space  freed  by 
compression  from  the  pressure  to  the  pressure  pc  and  the 
space  filled  by  the  incoming  liquid. 

When  a front  travels  to  the  right,  the  initial  free 
space  in  a layer  of  liquid  dx  just  back  of  the  front  is 
given  by 

Thereafter  we  have  a change  in  (j>  at  a rate  of  • 

Hence  at  any  time  t subsequent  to  the  passage  of  the  front 
at  the  time  t^,  we  have 

4>  s £ (^/3f ) t 

This  shows  that  the  value  of  (p  may  vary  from  point  to 
point  in  the  cavitated  region,  and  even  though  (p6  is  small, 
the  value  of  <p  may  increase  to  any  magnitude. 

4.  Motion  of  Plane  Breaking  Fronts 

No  influence  upon  the  molecules  of  a liquid  can  be 
propagated  with  a velocity  greater  than  the  velocity  of  sound 
in  the  liquid  past  a breaking  front.  Hence  the  pressure  and 
particle  velocity  ahead  of  a front  is  determined  by  the 
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initial  and  boundary  conditions  elsewhere.  If  the  boundary 
of  the  cavitated  region  moves  with  a velocity  less  than  the 
speed  of  sound,  then  cavitation  is  able  to  influence  the 
pressure  and  particle  velocity  in  its  neighborhood  where  the 
liquid  is  uncavitated.  Any  one  dimensional  disturbance  in 
an  unbroken  liquid  can  be  resolved  into  trains  of  waves 
going  in  opposite  directions  at  the  same  speed.  If  we  let 
Pi,  p2  and M(,M2  denote  corresponding  pressures  and  particle 
velocities  respectively,  then  we  can  write  for  the  two  waves 


Pi  = f C M, 

P - Pi  + P2 
Pi  - i(p  ifc .ajl) 


P2  - - 9 C M i 

^ i t *AA  j 

P2  " s(p 


If  we  substitute  these  into  the  condition 


> 0 


then 


> 0 


If  substituted  into  the  condition 


and  both  sides  are  squared,  then 
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The  first  and  last  of  the  preceding  inequalities  can  be 


true  if^^'J'^is  positive  and is  negative.  Hence  the 

conditions  for  the  propagation  of  a plane  breaking  front, 
along  with  the  expressions  for  and  mc  , can  be  written, 
since  in  the  form 


From  the  preceding  conditions  we  can  conclude  that 
any  wave  of  tension,  such  as  an  ultrasonic  wave,  can  give 
rise  to  cavitation  only  when  it  is  superposed  on  a wave  of 
decreasing  pressure  which  is  travelling  in  the  opposite 
direction.  This  is  exactly  what  happens  when  an  ultrasonic 
stationary  wave  is  set  up  in  a liquid. 

5.  Plane  Resting  Boundaries  and  Closing-Fronts 

We  shall  now  consider  the  behavior  of  the  boundary  of 

4 

a cavitated  region  when  it  cannot  advance  asA plane  breaking 
front.  This  is  the  case  if  a liquid  is  cavitated  ultra - 
sonically.  Although  a cavitated  region  will  advance  as  a 
plane  breaking  front  in  a liquid,  it  is  a qualitative  fact 
that  the  motion  of  this  front  ceases  when  it  reaches  the 
nodes  of  a stationary  ultrasonic  wave  pattern.  The  boundary 
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may  stand  still  or  it  may  advance  into  a cavitated  region 
leaving  the  liquid  behind  it  uncavitated,  in  which  case  the 
boundary  is  called  a closing-f ront . Which  alternative  occurs 
depends  upon  the  magnitude  of 


p-'  * the  pressure  in  the  compressional 

wave  train  that  approaches  the  boundary 
in  the  uncavitated  liquid 

Me-  particle  velocity  on  the  cavitated  side 
of  the  front  boundary 

the  gradient  of  the  particle  velocity 

e Y' 

( p - the  fraction  of  the  space  occupied  by 
the  cavities 


In  general,  we  have  three  cases  of  importance 
(a)  Intrinsic  closing-fronts 

In  this  case  we  have  for  initial  conditions 


where  x is  the  distance  from  the  boundary  into  the  cavitated 


Because  of  the  first  two  conditions,  cavitation  dis- 
appears in  the  layers  next  to  the  boundary;  hence  the  bound- 
ary advances  into  the  cavitated  region  as  a closing-f ront, 
leaving  the  uncavitated  liquid  behind.  This  is  called  an 
intrinsic  closing-front  and  cannot  move  through  the  liquid 
at  velocities  greater  than  c. 


and 


region. 
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(b)  Resting  boundaries 

Here  we  have  the  condition 

- ft  + S c 

The  pressure  is  continuous  at  the  boundary  and  equal  to  pc. 
When  the  above  equation  holds  but  the  boundary  cannot  ad- 
vance as  an  intrinsic  closing  front  at  a velocity  greater 
than  c,  then  the  boundary  remains  stationary  and  the  inci- 
dent waves  reflected  from  it  are  described  by  the  condition 
that 


where  p"  is  equal  to  the  pressure  in  the  receding  wave  train. 

(c)  Forced  closing-fronts 

The  initial  conditions  in  this  case  are 

, 0 > 0 

When  the  above  is  true  the  incident  waves  are  strong  enough 
so  that  the  boundary  advances  into  the  cavitated  region  and 
causes  the  bubbles  to  collapse.  This  action  causes  a finite 
discontinuity  of  pressure  at  the  boundary  if  4>>0.  A 
boundary  moving  in  this  fashion  is  called  a forced  closing- 
front.  It  can  be  shown  that  these  travel  at  a speed  less 
than  c relative  either  to  the  unbroken  liquid  behind  them  or 
to  the  cavitated  liquid  in  front. 
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6.  Three  Dimensional  Cavitation 

By  generalizing  from  the  one  dimensional  equations 
stated  above,  we  get  for  the  acoustical  equations  describing 
three  dimensional  phenomena 

if  ; -ff  - 

whe  re 


i,  j,  k are  unit  vectors  in  the  x,  y,  z directions 
respectively,  and 

= 2*3 L + 2±!Jj-  ■+■  ih!x 

^3 

"Cavitation  will  begin  at  a point  where  a local  minimum  of 
pressure  sinks  to  the  breaking  pressure  p^,  and  from  this 
point  a breaking-front  will  spread  out  in  the  form  of  a 
closed  surface  surrounding  the  point  of  initiation  and  en- 
closing a cavitated  region."  (33) 

Generalizations  from  the  one  dimensional  treatment 
lead  to  conditions  to  be  satisfied  in  the  uncavitated  liquid 
so  that  a front  may  advance,  expressable  as 

faiUrj  > 0 - &)  < fc 

where  •/K  is  equal  to  the  distance  into  the  unbroken  liquid 
along  the  normal  to  the  boundary.  Also  for  U^,  the  velocity 


(33)  E.H.Kennard,  p.180,  Ref. 32 
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of  a front  along  its  normal  relative  to  the  unbroken  liquid 
ahead,  we  get 


> C 


By  similar  generalizations  we  get  for 


In  general,  great  difficulty  is  encountered  when  a complete 
theory  of  three  dimensional  cavitation  is  attempted  hydro- 
dynamically . Consider  for  example  the  special  case  of 
spherical  waves. 

Let  P]_  f (ct  **  r)  be  the  pressure  in  a train  of 

waves  from  a fixed  center,  where  r is  a parameter  from  that 
center.  Then  let  these  waves  be  superposed  upon  other  waves 
expressed  as  P2(x,y,z,t).  In  this  case 


o 6. 


iP  C-  <Q/Cv  sMs 


can  be  written  by  means  of 


= ~ ?C  * JLat  aa* 


and  p * P]_  «#-  p2  as 

— Api — ^pT,  ^ d [Apj-  0 

If  the  normal  makes  an  angle  © with  the  radius,  then 

^ a = ~~  (*  + -£t)  *++  & 

Hence  the  necessary  condition  for  the  propagation  of  a 
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breaking  front  is 

- ( I - Ox  &)  + C C<X&  > + C 

^ ^ ^-6  3 /vs 

From  the  above  it  can  be  shown  that  a breaking  front 
cannot  follow  a diverging  spherical  wave  of  tension  into 
quiescent  liquid  where  (p-,  <0,  ^ “ 0)>  f°r  then  we  have 
© -0)  since  over  the  front  p - - p^.  In  order  to  cause 

cavitation,  spherical  waves  of  tension  must  overrun  other 
waves  of  such  character  that  the  right  side  of  the  above 
equation  has  large  negative  values  much  greater  than  is 
actually  possible. 

7.  Thermodynamic  Considerations 

The  above  theory  is  based  solely  upon  the  laws  of 
conservation  of  matter  and  conservation  of  momentum.  The 
theory  reveals,  however,  that  if  Pk<pc,  then  the  passage  of 
either  a breaking  front  or  a forced  closing-front  is  very 
likely  to  leave  behind  it  a greater  amount  of  energy  than 
the  kinetic  and  ootential  energy  of  the  cavitated  liquid  can 
account  for.  It  is  very  possible,  as  has  been  suggested  by 
Smith’s  theory  of  pulsating  bubbles,  that  the  excess  energy 
is  dissipated  or  radiated  away  in  consequence  of  the  local 
oscillations  associated  with  the  impulsive  production  and 
destruction  of  bubbles. 

The  difference  between  p^  and  pc  implies  an  inherent 
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instability  of  liquids.  There  exists  no  hydrodynamic  or 
thermodynamic  reason  why  cavitation  cannot  occur  at  pressures 
lower  than  pc.  From  the  point  of  view  of  this  theory,  this 
leads  to  breaking  fronts  advancing  with  velocities  greater 
than  c, which  is  very  improbable  from  a region  already  cavi- 
tated.  If  we  attach  the  problem  of  cavitation  statistically, 
rather  than  hydrodynamically , we  arrive  at  relationships 
which  are  more  readily  verifiable  by  experiment* 
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STATISTICAL  THEORY  OF  CAVITATION 


1.  Introduction 

A liquid  standing  at  room  temperature  may  be  described 
physically  as  being  in  a certain  equilibrium  state.  If  the 
liquid  is  placed  in  an  ultrasonic  field,  the  equilibrium 
configuration  which  the  molecules  of  the  liquid  possesssed 
at  room  temperature  is  altered  until  a new  equilibrium  state 
is  reached.  It  can  be  seen  that  in  passing  fromone  state  to 
another  there  is  a definite  rearrangement  of  the  molecules 
of  the  liquid.  The  rate  at  which  this  rearrangement  pro- 
ceeds and  the  new  equilibrium  configuration  which  is  reached 
is  dependent  upon  the  temperature  of  the  liquid  and  a char- 
acteristic energy  of  the  molecules.  It  is  evident  that  an 
initial  equilibrium  configuration  mu3t  be  assumed.  From 
this  initial  assumption  a statistical  theory  of  cavitation 
can  be  developed. 


2.  Theory  of  Rate  Processes^^ 

The  Arrhenius  equation  showing  the  variation  of  the 
specific  rate  of  a physical  reaction  (such  as  a rearrangement 
of  the  molecules)  with  temperature  can  be  expressed  as 


k - Ae 


s/rt 


(34)  S.  G-lasstone,  K.J.Laidler,  H.  Eyring,  The  Theory  of 
Rate  Processes,  McG-raw  Hill,  New  York  (1941) 
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where  k - specific  rate  of  the  reaction 

A = collision  number  or  frequency  factor 
E = difference  in  heat  content  between  the 
activated  molecules  in  the  ultrasonic 
field  and  the  "inert"  molecules  at  room 
temperature 

The  quantity  E is  usually  termed  the  "energy  of  activation" 
of  a reaction.  This  represents  the  energy  that  a molecule 
in  the  "inert"  state  must  acquire  from  an  interaction  with 
the  ultrasonic  field  before  it  can  take  part  in  the  rearrange- 
ment to  the  final  configuration. 

It  has  been  suggested  by  Sohay^-^)  that  there  exists 
the  possibility  of  a dualism  in  the  realm  of  sound  energy 
analogous  to  the  dual  nature  of  electromagnetic  radiation. 

If  we  assume  the  basis  of  present  quantum  theory,  it  is 
improbable  that  a quantum  of  sound  energy  does  exist  in  the 
interaction  of  sound  and  matter.  This  is  borne  out  by  the 
fact  that  the  relatively  low  frequencies  attainable  by 
present  day  ultrasonic  generators  would  give  us  energies  far 
below  the  observable  activation  energies  necessary  for  a 
physical  reaction  to  take  place.  It  must  be  mentioned,  how- 
ever, that  this  concept  has  not  been  fully  explored. 

A physical  reaction  is,  as  has  been  pointed  out,  char- 
acterized by  an  initial  configuration  which  passes  over  by  a 
continuous  change  of  coordinates  into  a final  configuration. 
"There  is,  however,  always  some  intermediate  configuration 
which  is  critical  for  the  process,  in  the  sense  that  if 

(35)  B.K.Sohay,  Journ.  Acoust.  Soc.  Am.  17,  285  (1946) 
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this  system  is  attained  there  is  a high  probability  that 
the  reactions  will  continue  to  completion.  This  critical 
configuration  is  called  the  "activated  complex"  of  the  re- 
action, and  is,  in  general,  situated  at  the  highest  point  of 
the  most  favorable  reaction  path  on  the  potential  energy 
surface . We  may  regard  this  activated  complex  as  an 

ordinary  molecule  which  possesses  all  the  usual  thermo- 
dynamic nroperties,  except  that  motion  is  restricted  along 
the  reaction  coordinate  leading  to  decomposition  at  a def- 
inite rate.  It  is  now  possible  to  derive  the  concentration 
and  rate  at  which  these  activated  complexes  pass  through 
the  critical  configuration  of  the  activated  state  by  statis- 
tical methods.  The  product  of  the  concentration  and  rate  is 
equal  to  the  specific  reaction  rate.  The  important  concept 
to  bear  in  mind  at  this  point  is  that  there  exists  an  acti- 
vated complex  which  is  the  critical  intermediate  in  every 
rate  process. 

If  we  let  C*  represent  the  number  of  activated  com- 
plexes per  unit  volume  of  liquid  lying  in  a length  d repres- 
enting the  activated  state  at  the  top  of  a given  potential 
barrier,  and  v the  mean  velocity  of  crossing  this  potential 
barrier,  then  we  have  V/&  equal  to  the  frequency  at  which 
the  activated  complexes  surmount  the  potential  barrier. 

Hence  we  have  for  the  rate  of  reaction 

Rate  of  reaction  = G 1 _ v 


(3 6)  S.G-lasstone,  et  al,  p.10,  Ref. 34 
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We  can  treat  these  activated  Comdexes  as  ordinary  molecules 
if  we  write  their  concentration  at  the  top  of  a potential 
harrier  as 


h 

where  m*  is  the  effective  mass  of  the  activated  complex  in 


the  partition  function  for  translation  in  the  reaction  path 
(only  possible  motion);  this  partition  function  represents 
the  probability  of  the  occurrence  of  the  activated  complex 
at  the  top  of  the  potential  barrier.  By  classical  methods 
it  can  be  shown  that 


From  this  it  can  be  derived  that  if  k is  the  specific  reaction 
rate,  then 


where  K'  is  the  constant  for  the  equilibrium.  The  constant 
K*  may  be  related  to  AF'  the  standard  free  energy,  AH'  the 
heat  content  change,  and  A S’  the  entropy  change  by  the  usual 
thermodynamic  relationships.  Hence  we  may  now  write  the 
expression  for  the  specific  reaction  rate  for  activated 
complexes  surmounting  a potential  barrier  in  the  form 


the  coordinate  of  reaction.  The  factor  (2iT  m*  KT)2  d/h  is 


Rate  of  reaction  " C KT 

h 


k = KT  K 

h 


k = KT  e~  A F ' /RT 


h 


or 
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h 

It  will  now  be  advantageous  to  transform  the  equili- 
brium constant  K*  and  the  thermodynamic  functions  in  terms 
of  the  partition  functions  of  the  species  involved  by  the 
methods  of  statistical  mechanics.  If  this  is  done,  we  obtain 

k = KT.  F’  e'£  /RT 

h Tf. 

X- 1 

for  a system  of  N molecules  where  F1  is  the  partition  function 
of  the  activated  complex  per  unit  volume,  F^  is  the  partition 
function  of  the  i»  reactant  in  the  physical  rearrangement 
and  E is  the  difference  at  the  absolute  zero  between  the 
energy  per  mole  of  activated  complex  and  the  sum  of  the 
energies  of  the  reactants. 

We  have,  at  this  point,  formulated  enough  basic  theory 
to  explain  to  a degree  the  process  of  ultrasonic  cavitation 
in  liquids.  It  is  in  order  here  to  assume  that  the  reactants 
mentioned  above  are  the  molecules  of  a liquid;  the  potential 
barrier  which  is  to  be  surmounted  may  be  explained  in  terms 
of  the  cohesive  strength  of  the  liquid;  the  energy  put  into 
the  system  originates  in  the  ultrasonic  transducer;  the  initial 
configuration  is  the  "non-cavitated"  liquid;  and  the  final 
configuration  is  the  cavitated  liquid. 
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( 37) 

3.  Eyring's  Liquid  Structure  Theory 

Eyring  regards  a liquid  as  made  up  of  "holes"  nearly 
molecular  in  size  moving  about  in  matter.  These  holes  are 
considered  as  playing  the  same  part  in  a liquid  as  molecules 
do  in  a gas.  The  holes  are  assumed  to  be  a part  of  the 
natural  basic  structure  of  the  liquid.  This  concept  differs 
somewhat  from  Furth’s  theory  in  that  Fiirth  assumes  the 
natural  holes  in  a liquid  to  be  under  high  vapour  pressure 
from  within,  while  Eyring's  natural  cavities  are  devoid  of 
matter,  at  least  in  the  initial  equilibrium  state. 

Suppose  now  we  have  N molecules  which  make  up  a liquid. 
Each  of  these  N molecules  is  bound  to  the  others  by  bonds 
adding  up  to  the  total  energy  E which  is  given  by 

E = ‘S.  ^ €x 

where 

^ - number  of  bonds  of  the  1—  kind 
Ej_  r bond  strength  of  the  i^  bond 

To  vaporize  the  N molecules  requires  exactly  NE/2  energy, 
since  each  bond  belongs  to  two  molecules.  Providing  no  hole 
is  left  in  the  liquid  upon  vaporization,  it  then  takes  an 
amount  of  energy  equal  to  E/2  to  vaporize  one  molecule. 
However,  if  we  vaporize  one  molecule,  and  leave  a hole,  we 
must  supply  energy  equal  to  E.  If  we  return  this  molecule 
to  the  liquid  we  bring  back  energy  equal  to  E/2.  It  is  seen 
then  that  it  requires  the  same  energy  e/2  to  make  a hole  in 


(37)  H. Eyring,  Ref. 25 
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a liquid  of  molecular  size  as  it  does  to  vaporize  a molecule 
without  leaving  a hole. 

If  we  consider  a molecule  in  empty  space  it  is  evident 
that  it  has  a higher  entropy  than  when  in  a liquid,  since 
it  has  a large  number  of  accessible  states  without  increase 
in  energy.  The  same  is  true  of  the  hole  dissolved  in  the 
liquid.  Now  a hole  can  have  many  positions  in  a liquid  and 
still  have  the  same  energy.  These  holes,  like  molecules, 
can  be  pictured  as  having  a random  motion  with  a velocity 
characteristic  of  the  molecules  at  that  temperature.  A liquid 
can  be  thought  of  as  having  a definite  lattice  arrangement 
somewhat  like  a crystal.  Wherever  a hole  exists,  a molecule 
close  by  can  jump  into  an  empty  lattice  point,  which  leaves 
a hole  behind,  and  as  this  process  goes  on  continuously  each 
hole  contributes  a new  degree  of  translation  to  the  liquid. 

If  the  liquid  is  at  a temperature  below  the  critical  or 
boiling  temperature,  this  translation  in  the  liquid  is  not  as 
free  as  in  the  vapour  phase,  hence  the  partition  function 
for  holes  is  smaller  than  for  vapour  molecules.  Also,  below 
this  critical  temperature  a molecule  in  going  from  one  hole 
to  another  has  to  surmount  the  potential  barrier  mentioned 
above.  At  the  critical  temperature  the  mobilities  and  parti- 
tion functions  of  the  holes  and  the  molecules  are  equal,  hence 
for  holes  in  liquids  or  molecules  in  empty  space  we  have  the 
theorem  of  equal  entropies  and  energies.  It  is  intuitive 


' 


j . . ' v v t c >t  " OlkTOZL.  li  ‘LiiJLb  'V.O  ,-jY/  u 

■ f . r f ' •-1C)-  o ■.  - O ■ 

- • 

. 

■ fi  : V..  •,  ' ■ cf  : 

j:  o.  is  :j'  ordJ  ' ' ' v£rl  ’if.  d:v  r 

. 

r.L  ■ ■ : ■ ■ : 

. ' d 


ncl,  . 

, . ' *.'•?  ■ : 

■ 

■ J &X£fCJin  u 
« I 

d ■ • tup  l t i . . 

. 


! ‘ . - 


50. 


that  the  energy  required  for  the  association  or  dissociation 
of  holes  of  molecular  size  depends  primarily  upon  the  nature 
of  the  liquid. 

Since  the  physical  process  of  cavitation  takes  place  in 
both  the  liquid  and  gas  phase,  we  assume  for  a first  approxi- 
mation that  association  or  dissociation  does  not  affect  the 
distribution  ratio  of  holes  between  the  liquid  and  gas  phase. 
If  this  is  assumed,  we  can  say  that  at  the  critical  tempera- 
ture we  should  have  as  much  matter  in  a unit  volume  of  other- 
wise empty  space  as  we  have  holes  in  a unit  volume  of  matter. 
As  the  temperature  is  lowered,  a plot  of  the  average  density 
of  the  liquid  as  a function  of  temperature  would  almost 
parallel  the  temperature  axis,  the  density  increasing  slightly 
as  the  temperature  decreases.  Hence  the  average  density  of 
the  vapour  and  liquid  phase  should  be  nearly  independent  of 
the  vapour  pressure  or  the  temperature.  "Now  the  law  of 
rectilinear  diameter  states  that  for  actual  cases  such  a plot 
of  average  density  is  a straight  line,  but  that  it  may,  and 
in  general  is,  a little  inclined  to  the  temperature  axis. 

Thus  the  law  of  Cailletet  and  Mathias  follows  from  statistical 
mechanical  considerations  of  holes. 

The  theorem  of  equal  energy  and  nearly  equal  entropies 
for  holes  in  matter  as  for  molecules  in  empty  space  frequently 
applies  in  the  same  way  to  the  solid  state  as  to  the  liquid 
state.  We  can  surmise,  therefore,  that  the  vapour  is  a rough 


(38)  H.Eyring,  P.285,  Ref. 25 
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measurement  of  the  number  of  holes  in  the  condensed  phase. 

The  theory  of  holes  is  of  importance  in  all  phenomena  in 
which  a part  of  the  matter  in  a condensed  phase  moves  with 
respect  to  other  parts.  This  follows  from  the  fact  that  at 
the  borders  of  holes  such  motions  can  be  set  up  with  a mini- 
mum of  disturbance  to  molecules  not  participating  in  the 
motion. 

In  certain  phenomena  it  may  be  much  simpler  to  fix 
attention  on  the  few  holes  rather  than  on  the  much  more 
numerous  molecules.  Only  in  very  precisely  ordered  arrange- 
ments, can  the  attainment  of  a critical  motion  by  a few  mole- 
cules profoundly  affect  the  stability  of  neighbouring 
molecules  and  these  in  turn  release  still  others  from  their 
restraining  potential. 

Some  support  for  the  general  point  of  view  with  res- 
pect to  holes  as  developed  here  may  be  obtained  from  dis- 
cussions of  the  Wilson  theory  of  semiconductors.  In  this 
theory,  electrical  conductance  is  assumed  to  arise  from 
thermal  fluctuations  which  raises  an  electron  of  low  energy 
up  to  the  conduction  levels  leaving  a positive  hole  behind. 
This  theory  seems  to  account  satisfactorily  for  the  observed 
facts. < 39 > 

4.  Interpretation  of  Cavitation  and  Conclusions 

It  was  mentioned  earlier,  in  conjunction  with  Briggs’ 


(59)  A. Fowler,  Proc.Roy .Soc.  A140,  505  (1935) 
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experiments  on  the  cavitation  strength  of  liquids,  that 

the  only  property  of  a liquid  consistently  related  to  the  co- 

. ( 40 ) 

hesive  pressure  was  its  viscosity.  An  extension  of  Eyring  sv  ' 
theory  of  absolute  reaction  rates  to  viscosity,  coupled  with 
his  "theory  of  holes"  as  applied  to  liquid  structure,  yields 
some  information  concerning  this  relationship. 

If  we  consider  two  layers  of  molecules  a distance  S, 
apart,  with  a force  F0  tending  to  displace  one  layer  with 
respect  to  the  other,  then  the  coefficient  of  viscosity  is 

given  by  A --  AJi 

A <v 

where  v is  the  difference  in  the  velocity  of  the  two  layers. 

The  motion  of  a layer  of  molecules  is  assumed  to  involve 
passage  of  a molecule  from  one  equilibrium  position  to  another. 
For  this  passage  to  occur,  it  is  required  that  a suitable 
hole  be  available,  and  the  production  of  such  a hole  requires 
energy.  Hence  the  passage  of  a molecule  from  one  equili- 
brium position  to  another  is  the  same  as  the  system  surmount- 
ing a potential  barrier.  In  order  to  determine  the  energy 
involved,  it  is  convenient  to  adopt  the  following  notation. 

Let  S be  the  distance  between  two  equilibrium  positions  and 
the  distance  between  two  molecules,  both  measured  in  the 
direction  of  motion.  Also  let  be  the  mean  distance 
between  two  adjacent  molecules  in  a direction  at  right  angles 
to  the  direction  of  motion.  If  we  assume  a symmetrical 


(40)  S.  G-lasstone,  et  al., 
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potential  energy  barrier,  then  the  distance  between  the 
equilibrium  position  and  the  activated  state  (at  the  top  of 
the  barrier)  is  equal  to  s S • The  energy  acquired  by  a 
moving  molecule  when  it  has  reached  the  top  of  the  potential 
barrier  is  then  given  by  ■§•  F0  • 

With  the  assumption  that  the  transmission  coefficient 
for  the  passage  of  a molecule  over  a potential  energy  barrier 
will  be  unity  and  neglect  tunneling  through  the  barrier,  then 
the  rate  of  passage  of  a molecule  over  the  barrier  is 


while  the  reaction  rate  for  a molecule  that  has  once  sur- 
mounted a potential  barrier  to  return  to  its  initial  state 
may  be  stated  as 


k = KT  F'  e~E/RT 


h F 


( 41) 

when  no  force  is  applied.  According  to  Briggs  ' , the 
forward  reaction  rate  k is  given  by 


i ~ & S a S3  ) /&KT 

% e 

Then  A v for  each  molecule  in  the  faster  layer  is 


or 


A at  - £ S h.  ( Fo  S Sa  S3 


(41)  H.B. Briggs,  et  al,  p.674.  Ref.  13 
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With  this  equation  for^W  we  may  write  the  expression  for 
the  viscosity  as 

*1  - F»%,/  2%k  *^X(F0  Si*i3)/s/<T 

or  * * s.#T/sas4sa* 

if  we  assume  that  ZoSSaSa  ^ ^ 2KT,  as  is  the  case  for 
very  small  forces.  If  we  introduce  the  value  of  k and  Sj 
assume  $ equal  to  5,  , we  get  for  the  viscosity 

A Ue£/KT 

\J(l-  q-AiuMt) 

or 

4 = A e E/KT 

where  N is  Avogidro ' s number  and  V the  molal  volume  of  the 
liquid  which  may  be  taken  on  the  average  to  be  about  40  cm^. 
In  order  to  arrive  at  the  above  expression  for  the  viscosity, 
Eyring  assumes  that 

= (i  - e~i7</*T) 

where'll  are  frequencies,  having  a range  between  1.4  and 

1.8  x 10  , which  can  be  obtained  from  the  specific  heats 

of  the  liquids.  With  this  range  we  get  F'/F  values  between 

0.1  and  0.25.  From  the  average  value  of  V we  see  that  the 

factor  Nh/V  is  of  the  order  of  10”^.  This  must  be  multiplied 

by  a factor  4 to  10  to  account  for  f/F1.  Hence  the  constant 

i -4 

A in  the  last  equation  for  viscosity  varies  between  4 x 10 
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and  10"^,  depending  on  the  liquid.  We  must  now  consider 
the  activation  energy,  E,  which  is  related  to  the  internal 
energy  plus  the  energy  necessary  to  overcome  the  external 
pressure . 

Eyring  assumes  that  following  van  der  Waals'  equation 
(p-f-a/v2)  (v  - b)  - RT 
We  may  write  for  the  activation  energy 

E 3 ( p * a/v2 ) Q 

where  S is  a volume  equal  to  about  2v.  If  this  is  sub- 
stituted in  our  last  expression  for  the  viscosity,  we  get 

H*  4 e*p[(?  *R^y2)  ] 

When  the  pressure,  p,  is  zero  then 

E = (a/v2)  ( 5/N) 

It  is  clear  that  as  p increases  E increases  and  hence 
the  viscosity  should  increase.  If  we  have  a negative  pressure, 
which  is  always  the  case  when  cavitation  occurs,  the  energy 

r\ 

will  decrease,  and  if  p = -a/v  the  activation  energy  should 
decrease  to  zero  and  the  potential  barriers  disappear.  At 
this  point  the  liquid  becomes  gas-like;  p cannot  then  be 
more  negative  than  a/v2,  the  van  der  Waals  cohesive  pressure. 

A 

As  was  mentioned  earlier,  a negative  pressure  of  about  10 
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atmospheres  fulfills  this  limiting  condition.  This  is  about 
of  the  same  order  of  magnitude  as  the  theoretical  cohesive 
pressure  for  a liquid. 

However,  by  similar  considerations  for  solids  the  actual 
strength  has  been  found  experimentally  to  be  only  1/1000  of 
the  theoretical  values.  Then  for  liquids  we  should  expect 
to  have  to  multiply  p by  a factor which  is  of  the  order  of 
1000  when  the  applied  pressure  is  negative.  If,  however,  the 
pressure  is  positive,  then  it  is  very  probable  that  the  factor 
is  close  to  unity  in  order  to  get  the  observed  dependence 
of  viscosity  and  pressure.  Since  in  an  ultrasonic  field  the 
pressure  varies  with  time  from  positive  to  negative  values, 
we  should  expect  /3  to  be  a function  of  time.  Then  we  can 
wri  te 

fi  - A e*P<  i/i  M aA2?/i?  ) 

1 ( RT  ) 

where  ^ is  a function  of  time  and  the  ratio  p/pc  (pc  being  the 
cohesive  pressure  of  the  liquid)  is  equal  to  1000  for  steady 
state  negative  pressure  nearly  equal  to  the  cohesive  pressure, 
and  also  equal  to  unity  for  positive  pressures. 

If  we  assume  the  theoretical  derivation  of  these  equations, 

then  we  can  interpret  some  of  the  experimental  results  obtained 
(42) 

by  Briggs.  In  Fig.  4 we  have  plotted  the  cavitation  for 

degassed  kerosene.  The  value  of  a/v  for  this  liquid  is  small 


(42)  H.B. Briggs,  et  al,  p.675>  Ref.  13 
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but  positive,  and  it  was  found  that  it  requires  a negative 
pressure  of  about  0.8  to  1 atmosphere  to  cause  cavitation. 
With  the  application  of  a hydrostatic  pressure  we  change 
a/v  but  little,  and  hence  it  requires  the  same  negative 
pressure  to  cause  cavitation.  However,  since  we  have  intro- 
duced a hydrostatic  pressure  of  A atmospheres,  it  requires 
an  acoustic  pressure  of  A + 1 atmospheres  to  produce  a 
negative  pressure  in  the  liquid  of  one  atmosphere.  We 
obtain,  therefore,  a straight  line  relationship  between  the 
amplitude  of  the  sound  pressure  and  the  hydrostatic  pressure 
of  the  liquid.  We  assume  here  that  cavitation  input  voltage 
is  a function  of  the  sound  pressure.  The  straight  line  ob- 
tained is  displaced  from  the  zero  hydrostatic  pressure  origin 
by  the  amount  of  the  cohesive  pressure  of  the  liquid. 

The  last  equation  for  the  viscosity  of  a liquid,  derived 
above,  explains  the  dependence  of  cavitation  pressure  on  vis- 
cosity observed  in  Fig. 5 and  given  by  the  equation 


3.osfc 


If  we  assume  that 


A = 


Nk , , 

v (i  - e- Wat] 


is  equivalent  to  ^ above,  then 

Ju(yA)  = AS/rTv'  --  E/ XT 

where  is  fairly  constant  from  liquid  to  liquid. 

Also  if  we  assume  that  the  value  of  $ is  constant  from 
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liquid  to  liquid,  we  should  have  cavitation  occurring  at  a 
negative  pressure  equal  to 


“Pc 


Hence  we  may  now  evaluate  the  constant  since 


or 


3.05  pc 

10  6 

^Ja  * 

(3.05  a/v2) 

( 106  /3  ) 

Va- 

(a/v2)  S 
RT 

/3  - 

(3.05)  RT 
io6  S 

if  pc  is  given  in  dynes/ cm2.  Now  RT  at  300°K  is  equal  to 
2.42  x lO"1"0,  and  according  to  Eyring  S is  about  2v  = 80  cm-'1. 
Hence  is  about  930,  which  is  in  good  agreement  with  experi 
mental  observations. 

The  curve  of  Fig.  & indicates  that  the  constant  /£  is  a 
function  of  time.  If  we  take 

/3  = f 30  ( I - 3)  + i 

where  o(  is  about  2.2  for  castor  oil,  then  the  experimental 
curve  is  fitted  within  the  accuracy  of  measurement.  It  can  be 
assumed  that  the  above  equation  holds  down  to  very  short  time 
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intervals;  then  the  time  necessary  for  /j  to  double  is  given 

t>y 


to*  1/(930)  (2. 2) 
or 

t0  - 1.2  x 10  10sec. 


which  is  the  time  required  for  a force  of  half  the  theoreti- 
cal breaking  strength  to  produce  a cavity  in  the  liquid. 

"This  should  probably  occur  when  the  cavity  size  has  grown 
from  1 to  2 units,  and  hence  should  compare  with  the  trans- 

( 43) 

ition  time  for  a molecule  jumping  from  one  hole  to  another." 
Eyring  has  calculated  this  time  to  be 


t « 1_  = ( _£_)  ( F ) 
k ( KT)  ( W) 


eS/KT 


or 

t - 1.6  x 10”12  eE//KT  sec 


For  castor  oil  upon  which  no  external  force  is  acting,  we 
get  from  activation  energy  measurements  a value  for  E/KT 
equal  to  8.6.  Now  for  a negative  pressure  of  half  the  break- 
ing strength  this  factor  is  reduced  to  4.3  and  hence  the  time 
should  be 

tQ  = 1.6  x 10~12e"4*5  = 1.15  x 10“10  sec. 
which  is  in  good  agreement  with  the  value  of  1.2  x 10  E°sec. 
found  experimentally. 

(43)  H.B. Briggs,  et  al.  Ref.  42 
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It  would  seem  that  the  more  recent  experimental  evidence 
is  in  favor  of  Eyring's  statistical  theory  of  cavitation,  and 
that  a thorough  hydrodynamic  theory  becomes  too  complex  to 
be  interpreted  and  tested  experimentally. 
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ABSTRACT 

The  input  of  ultrasonic  power  into  a liquid  medium 
causes  a definite  rearrangement  of  the  initial  molecular 
configuration.  When  cavitation  occurs,  it  is  ooserved  that 
small  cavities,  filled  with  the  vapour  of  the  liquid,  are 
formed  in  the  liquid.  These  cavities  vanish  as  soon  as  the 
ultrasonic  power  source  is  shut  off,  and  the  liquid  presum- 
ably returns  to  its  original  state.  It  is  necessary,  there- 
fore, to  assume  an  initial  molecular  configuration,  i.e., 
postulate  a basic  liquid  structure  theory,  the  liquid  being 
under  standard  conditions  of  pressure  and  temperature. 

Further,  a "transition  theory"  describing  the  passage  of  the 
initial  molecular  configuration  to  the  final  cavitated  con- 
figuration must  be  worked  out.  Cavitation  in  a gassy  and  a 
completely  degassed  liquid  would  then  seem  to  require  two 
separate  theories  for  its  occurrence,  since  their  initial 
molecular  configurations  would  appear  to  differ  widely.  This 
difficulty  can  be  surmounted,  however,  if  we  assume  that  a 
completely  degassed  liquid  is  made  up  of  natural  holes  which 
are  molecular  in  size  and  even,  perhaps,  the  size  of  two  or 
three  molecules.  Holes  larger  than  the  size  of  a molecule 
may  be  formed  by  the  combination  of  several  unit  size  holes. 
This  produces  a metastable  condition,  unless  the  liquid  happens 
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to  be  in  a force  field  which  would  activate  and  sustain 
such  a chain  of  events.  Hence,  an  ultrasonic  field  gives  us 
a source  for  the  required  "activation  energy",  as  well  as  a 
source  of  "sustaining  energy"  resulting  in  further  combina- 
tion of  the  "nuclei"  cavities  until  a definite  break  or 
fracture  occurs.  Although  this  point  of  view  can  be  applied 
strictly  to  degassed  liquids,  it  is  probable  that  a completely 
degassed  liquid  does  not  exist.  Hence,  it  would  appear  more 
tenable  to  assume  that  both  natural  holes,  as  well  as  micro- 
scopic gas  bubbles,  exist  in  liquids  under  standard  conditions. 
Although  at  the  initial  cavitation  point  the  natural  holes 
and  the  gas  bubbles  behave  differently,  it  is  probably  a very 
short  time  before  the  cavities  which  cause  the  fracture  of 
the  liquid  consist  of  gas  bubbles  alone.  These  are  the  visible 
cavities  observed  in  cavitation  experiments. 

It  has  been  observed  that  cavitation  takes  place  during 
the  negative  half  of  the  pressure  cycle  of  an  ultrasonic 
field.  If  we  assume  that  a certain  activation  energy  must  be 
supplied  to  a molecule  to  cause  it  to  alter  the  initial  con- 
figuration, then  it  is  reasonable  to  conclude  that  the  effect 
of  a negative  pressure  is  to  decrease  the  potential  energy 
barrier,  making  it  easier  for  a molecule  to"jum]3'  away  from  a 
hole.  This  causes  the  cavities  to  increase  in  size.  The 
probability  that  the  effect  of  a negative  pressure  is  to  in- 
crease the  tunneling  through  the  potential  can  be  taken  as 
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being  too  low  to  have  any  effect.  Conversely,  the  effect  of 
a positive  pressure  is  to  increase  the  potential  carrier 
for  a molecule  to  jump  away  from  a hole,  tending  to  decrease 
the  size  of  the  existing  cavities. 

If  the  above  reasoning  is  followed,  it  appears  that  a 
reversible  process  is  occurring,  i.e.,  we  have  the  formation 
of  cavities  during  the  negative  half  of  the  pressure  cycle, 
and  the  destruction  of  cavities  during  the  positive  half  of 
the  cycle.  This  must  be  true  so  long  as  the  sound  pressure 
does  not  exceed  the  cohesive  pressure  of  the  liquid.  How- 
ever, once  the  cohesive  pressure  is  exceeded  the  reversibility 
ceases,  with  the  result  that  cavities  are  formed  and  increase 
in  size  until  there  is  complete  rupture  of  the  liquid.  There 
is  also  the  possibility  that  the  existence  of  microscopic 
bubbles  in  the  liquid  may  cause  the  reversibility  to  shift 
from  equilibrium  at  some  point  before  the  sound  pressure 
exceeds  the  cohesive  pressure  of  the  liquid.  Recent  experi- 
ments with  liquids  known  to  contain  these  gas  bubbles  seem  to 
indicate  that  this  is  the  case. 

The  relationship  between  the  activation  energy  and  the 
cohesive  strength  of  a liquid  has  been  determined.  From  this, 
it  can  be  concluded  that  the  statistical  theory  of  cavitation, 
as  suggested  by  Eyring's  theory  of  physical  processes  in 
liquids,  is  the  more  tenable  theory  of  cavitation,  and  offers 
the  better  avenue  for  further  investigation. 
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Figure  2 
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Figure  4 
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Figure  5 
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Figure  6 
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TABLE  I 


Name 

Vise . 
in 

poises 

V.P.* 

'tr 

/ 

X 10 

? 

Pc 

Castor  Oil  D*B. 

6.30 

A 

1.477 

0.969 

3.90 

Olive  Oil 

0*84 

B 

1.431 

0.912 

3.61 

Corn  Oil 

0.63 

B 

1.463 

0.914 

3.05 

Peanut  Oil 

— 

B 

1*458 

0.936 

3.08 

Soya  Bean  Oil 

0.45 

B 

1*461 

0.919 

2.45 

Cottonseed  Oil 

— 

B 

1.463 

0.912 

2.32 

Linseed  Oil 

0.38 

B 

1.468 

0.921 

2.36 

Sperm  Oil 

0.25 

B 

1.440 

0.880 

2.90 

Kerosene 

0.04 

B 

1.324 

0.810 

2.00 

G-*  E*  Transil 

0.13 

B 

1*350 

0.880 

2.35 

* Vapor  pressure  ranges  A 0,01  cm*  Hg*,  B 0*1  cm*  Hg* 
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